THE ELECTROMAGNETIC MASS DIFFERENCE OF PIONS 

AT LOW TEMPERATURE 



Cristina Manuel 
Dpt. Estructura i Constituents de la Materia 
Facultat de Fisica, Universitat de Barcelona 
Diagonal 647, 08028 Barcelona (SPAIN) 

Q\ ; Nuria Rius 

. Dpto. de Fisica Teorica and IFIC 

Centro Mixto Universidad de Valencia-CSIC 
i 46100 Burjasot, Valencia (SPAIN) 



X 



Abstract 

We compute low temperature corrections to the electromagnetic mass dif- 
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independent way in the framework of chiral perturbation theory, using the 

background field method and the hard thermal loop approximation. We also 
S-H. 

D . 

generalize at low temperature the sum rule of Das et al. We find that the 



mass difference between the charged and neutral pions decreases at low tem- 
perature T with respect to the T = value. This is so in spite of the fact 
that charged particles always get a thermal correction to their masses of order 
~ eT, where e is the gauge coupling constant. Our result can be understood 
as a consequence of the tendency towards chiral symmetry restoration at finite 
temperature. 
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I. INTRODUCTION 



At low energies the strong interactions are successfully described in the framework of 
chiral perturbation theory (%PT) [|l],|2|]. This theory only involves the low energy modes of the 
QCD spectrum, such as the (pseudo) Goldstone bosons of the spontaneously broken chiral 
symmetry. There are eight (pseudo) Goldstone bosons, (it's, K's, 77), whose interactions can 
be understood in terms of symmetry considerations. A chiral Lagrangian is expanded in 
derivatives of the Goldstone fields, and in the masses of the three light quarks, which break 
explicitly chiral symmetry. Electromagnetic interactions can also be included in a chiral 
Lagrangian, although they break the chiral symmetry explicitly. 

The mass difference between the charged and neutral pions in the chiral limit (i.e. 
massless quarks) was computed more than thirty years ago using a sum rule approach 
||. With very general hypothesis, namely, soft pion theorems, current algebra, Weinberg's 
sum rules [Q , and vector meson dominance, Das et al. reached to the theoretical prediction 
(M n ± — M n o)th ~ 5.0 MeV. That value was in a remarkable good agreement with the ex- 
perimental one (M n ± — M v n) exp ~ 4.6 MeV. This computation allowed one to understand 
this mass difference as being due, essentially, to electromagnetic effects. This result was also 
reproduced in the framework of x?T, by considering resonance exchange within a photon 
loop ||. Better theoretical estimates of the electromagnetic mass difference of pions have 
since then been achieved 0. 

In the chiral limit Dashen's theorem [|7| states that the mass difference between charged 
and neutral kaons should be the same as the one of pions. Electromagnetic effects can be 
parametrized in a chiral Lagrangian in such a way that Dashen's theorem is fulfilled in the 
chiral limit |J. However, there is no good agreement between the predicted electromagnetic 
mass difference and the experimental value (M^±— M^-o) ea . p ~ —4.0 MeV. This is understood 
from the fact that the kaon's mass difference receives important QCD contributions of the 
order of the quark masses. On the other hand, corrections to Dashen's theorem due to the 
effect of the quark masses have been studied thoroughly in the literature (see, e.g. Refs. 
and || and references therein). 

In this article we will compute the electromagnetic mass difference of pions and kaons 
in the chiral limit at low T. Exactly as it happens at T = 0, we expect that our results 
for those masses should be trusted for pions. To get a better estimate of [M K ± — M K o){T) 
further corrections should be taken into account, but we will not compute them here. We 
will only consider low temperatures, that is, T -C f n , where f n = 92.4 MeV is the pion decay 
constant. This is the temperature regime where the contribution of the lightest particles 
of the QCD spectrum is the relevant one, since those of heavier states are exponentially 
suppressed. Thus, the thermal contribution of vector resonances will not be taken into 
account in this article. 

The computation will be first done in a model independent way in the framework of xPT, 
using the imaginary time formalism. We will consider the lowest order chiral Lagrangian 
extended to include electromagnetism, and obtain the one-loop thermal effective action for 
soft modes with the help of the background field method. We will work in the hard thermal 
loop (HTL) approximation [HI 111]. This methodic approximation allows us to identify the 
leading thermal corrections to any Feynman diagram, while ensuring the gauge independence 
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of the results. 

Finally, we will generalize to low T the sum rule computation of Das et al. ||. Using the 
same kind of soft pion techniques as in the T = case, we will check our previous results for 
pions. In this approach we will make use of the evaluation of thermal correlators of vector 
and axial vector currents derived in Ref. ||12|| . 

We will make a big emphasis on the assumptions and approximations involved in any of 
the two above mentioned approaches. It will then be clearly established that the standard 
soft pion techniques at low T in the chiral limit have their diagrammatic counterpart in the 
HTL approximation. 

This paper is structured as follows. In Subsect. II A we review the \PT framework used 
in our computation. In Subsect. [II B| we obtain the one-loop thermal effective action for the 
soft Goldstone bosons. The sum rule computation of Das et al @ is reviewed in Subsect. 
Ill A| . This computation is then generalized at low T in Subsec. [IIIB| . In Subsect. |III C 



we present the numerical results, and we end with the conclusions in Sect. IV. A general 
discussion on HTL's is presented in the Appendix. 



II. CHIRAL PERTURBATION THEORY COMPUTATION 



A. Lowest Order Chiral Lagrangian with Electromagnetic Interactions 

In this subsection we present a brief summary of in the lowest order extended to 
include electromagnetic interactions. We also explain how the background field method 
(BFM) is used in this specific case. We use natural units, so that Ti = c = = 1. 

At leading order the effective low energy Lagrangian describing the physics of the 
(pseudo) Goldstone bosons of the spontaneously broken chiral symmetry is, in Minkowski 
space-time B 



f 2 / N f 2 



(2.1) 



+ J -f Tr (v^VE) + ^Tr ( X f S + x£ f ) + CTr (QEQE 1 



where is the photon field, F^ v is the photon field strength, a is the gauge- fixing parameter, 
and E is a SU (N) unitary matrix written in terms of the (pseudo) Goldstone bosons $ as 



E = exp (i$/U) 



(2.2) 



For N = 3, = $ a A a , where A a are the Gell-Mann matrices, and in terms of the physical 
particles, one can write 



$ = V2 



/ 7T U I 778 



71 

\/2 + y/E 

K° 



K + \ 
K° 

_2?78 / 
V6 I 



(2.3) 



2 



The matrix Q is the quark charge matrix, which for the N = 3 case reads 



/2 



(2.4) 



V 



-1/ 



and e is the electromagnetic coupling constant. 
The covariant derivative is defined as 



QA 



(2.5) 



ty, and a M being external vector and axial vector hermitian and traceless sources, respectively. 
One also defines 



X = 2B(s + ip) 



(2.6) 



where B is a constant related to the quark condensate (qq) = —f 2 B[l + 0(m q )], and s 
and p are scalar and pseudoscalar external sources, respectively. The external scalar source 
can incorporate the mass matrix of quarks, which we put to zero since we will work in the 
chiral limit. In order to mantain the usual chiral counting, it is convenient to count the 
photon field as a quantity of order 0(1) and the electromagnetic coupling e of 0(p), so the 
Lagrangian £ 2 is of order 0(p 2 ) ||. 

The parameter C gives the purely electromagnetic part of the masses of the charged 
pions and kaons in the chiral limit 



M 2 



K± 



2e 2 C 
j - 



0(m q 



(2.7) 



The constant C obeys a sum rule which relates its value to the correlators of vector and 
axial vector currents (see Subsec. [Ill A[ ). The low energy constant C is a model independent 
quantity, which enters naturally in the chiral Lagrangian. 

To ensure the SU(N)r x SU(N)l symmetry of the Lagrangian (|2.1|) one introduces local 
spurions Q R (x) and Q L (x) instead of the constant charge matrix, Eq. fl2.4|) . The Lagrangian 
fl2.ll) is then modified so that the covariant derivative becomes 



(2.8) 



where, for convenience, we have defined 



v u + a 



A' ■ 



The last term in Eq. (|2.1|) also becomes 

Under a SU(N)r x SU(N)l symmetry all fields transform as 



(2.9) 



(2.10) 
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= Un(x)ll(x)U\ (x) 






(2.11a) 








= U l {x)Q l {x)U}{x) , / = 


R, L , 
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) + ip(x))' 
U t (x) 


= U R {x) (s(x) + [/|(^ 
e 5C/(iV)j , l = R,L . 


) , 




(2. lie) 
(2.11f) 



At this point one is ready to compute one-loop corrections to the lowest order chiral La- 
grangian with the background field method. Let us recall that at the end of the computation 
one has to put the spurions Qr(x) = Ql( x ) = Q- 

In the BFM all fields are split into classical and quantum pieces. The photon gauge field 
is split additively 

A^x) = A^x) + a^x) , (2.12) 

where A^ is the background field, and is the quantum one. The matrix containing the 
Goldstone fields is split multiplicatively 

where 

£(x) = , (2.14) 

is the background or classical field, and is the quantum fluctuation. 

In the spirit of the BFM, one expands the Lagrangian keeping only terms which are 
quadratic in the fluctuations 

£ 2 = 4 0) + 4 2) + --- (2.15) 

We will work in Euclidean space-time, so that we rotate the Minkowski Lagrangian £2 
to the Euclidean one Li,e- In Euclidean space-time one has |§ 

4=4^-^w a (2 - i6) 

- fl Tr(A M ) 2 - 4 Tr(M + ) - % Tr (H 2 R - R\ 



and 



4> = - i Tr ([A„0]j + -Tr(M+0 2 ) (2.17) 



4 - d,^) 2 - ^(d^) 2 + yTr ([A M , 0] iZ*) a" 

^Tr([/^ + /i L ,0][/^-tf L ,0]) 

^Tr (d„ a" + S-Tr{H 2 L )a^ 
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where 



d^4> = + [f „, 0] , (2.18a) 

r M = \ (e f vje + evje f ) , (2.18b) 

4. = ^ (e f - ev^ f ) , (2.1*0 

V\ = d^-i{Fl + Q l A^) , l = R,L, (2.18d) 

M + = ^ t + ^, (2-18e) 

fij i = e t ^+^ t , (2.isf) 

H l = ?Qr£-ZQl? , (2.18g) 



We will work in the massless quark limit. In this situation all particles at tree level are 
massless except for the charged pions and kaons. 

To obtain the one-loop effective action one has to integrate out the quantum fluctuations 

(f) and a M , while treating the background fields as external sources. The functional integral 

(2) 

can be done since C 2 E is only quadratic in the quantum fields. In a diagrammatic approach, 
one should only consider diagrams with quantum fields running inside the loop, while the 
background or classical fields only generate external vertices. At T = the one-loop func- 
tional generated by 4b has been evaluated in Ref. |§, using dimensional regularization 
and renormalization to deal with the ultraviolet (UV) divergencies. In Ref. || the next-to- 
leading order chiral Lagrangian £4 has also been constructed. We will not need it at the 
order of the computation we are working. 



B. One-Loop Thermal Effective Action for Soft Modes 

In this subsection we use the BFM to compute one-loop thermal corrections to the low 
energy constant C for soft modes. We will only consider thermal corrections and refer to || 
for a more general T = analysis. Let us mention that the UV divergencies which appear 
at T = and at finite T are the same, so the renormalization procedure does not change at 
all in our analysis with respect to the T = one. 

With respect to UV finite corrections in the chiral limit, and if we restrict our computa- 
tion to soft modes, that is momenta P <T, then one can ensure that the leading thermal 
corrections, which are typically of order T 2 /f 2 , are dominant with respect to the T = 
ones, which are of order P 2 /f 2 . Therefore, if one neglects corrections of order P 2 /f 2 , there 
are only two expansion parameters in the chiral limit. One of them is the gauge coupling 
constant squared, e 2 , and the other is the dimensionless quantity T 2 / f 2 . We will assume 
that the two expansion parameters are equally important. 

The computations will be done using the imaginary time formalism (ITF). Feynman 
rules for propagators and vertices are straightforwardly derived from Eq. (|2.17|) . We will 
denote Euclidean momentum with capital letters, so K 2 = k 2 , + k 2 . In ITF and for bosonic 
fields ko = 27rnT for integral n. The standard notation for the thermal momentum measure 
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lAjy- n=oo . J3t, 

= T H / ?t4 > ( 2 - 19 ) 



(2tt) 4 n ^J (27T 

will be used throughout this article. 

The leading thermal corrections to any Feynman diagram when the external momenta is 
soft, that is <C T, arises when the momenta running inside the loop is hard, or ~ T. Those 
diagrams are called hard thermal loops (HTL's) [pT0|,p7Tf| . There are systematic set of rules 
to extract from each Feynman diagram the corresponding HTL. Those will be taken into 
account in the computation. 

In this paper we are interested in finding the thermal corrections to the masses of the 
soft charged pions and kaons. So we are looking for one-loop corrections to the term in C^ E 

- jTr (H 2 R -H 2 ) . (2.20) 

Since we will only compute the leading thermal correction to the effective action, and 
the masses of the charged quantum fields are soft, we will neglect the masses of the quantum 
fields in the propagators. This means that we are neglecting corrections of order M n ±/T 
in the final answers. On the other hand, considering those masses as soft quantities means 
than we are assuming T 2 > Cj f 2 . 

There are only three Feynman diagrams to be evaluated, namely those generated by 
the last three terms of Eq. (|2.17|) , all of which give, in the HTL approximation, a tadpole 
type contribution at one-loop order. We use dimensional regularization to deal with the UV 
divergencies. In dimensional regularization and renormalization the T = contribution to 
the tadpole vanishes, so we only need to consider the thermal contributions to each diagram. 

The first diagram we consider is a tadpole of quantum fields 0. The generators of SU(N) 
\ a are normalized as Tr(A a A b ) = 2S ab , and [A a , \ b ] = 2if abc \ c , where f abc are the structure 
constants of SU(N). Then this tadpole gives the correction to Eq. ( |2.20|) 

C_ f a b c f a> bc{HR + HL) a {HR _ J^^=C^Tr(H 2 R - H 2 ) . (2.21) 



In Eq. ( P^TD we have used the SU(N) relation f abc f a ' bc = N6 aa ' 



Taking into account the relation 

TrH 2 L = -i Tr (h% - H 2 ) + Tr (Q 2 R + Q 2 L ) , (2.22) 



one sees that there are two other diagrams which give corrections to Eq. ( [2.20|) . One of 
them is a tadpole of the quantum gauge field, and the other is a diagram with a quantum 
boson (ft and a quantum gauge field circulating inside the loop. We evaluate the two 
above diagrams in the Feynman gauge a = 1. However, since we only retain the HTL of the 
diagram, one can ensure that our results are gauge independent [p~0| , p~T]| . We thus have 

(2.23) 




d A K 5^5^ r d 4 K (K + P)„(K + P) v 6^ * 



(2tt) 4 K 2 J (2tt) 4 K 2 (K + P) 2 
Using the relation (|2.22|) we end up with 
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-^Tr(H 2 R -Hl) —— = -^Tr(H%-Hl) , (2.24) 



flr^ /rr2 tt*\ f ^ K 3 f*T 2 

(2vr) 4 K 2 ~ 32 

plus an irrelevant constant term proportional to Tr (Q 2 R + Q 2 L ). 

The one-loop thermal effects for soft background fields involving the T M and A M fields 
arising from the two first pieces of Eq. ( 2.17]) have already been obtained in Refs. flT3| , |T4" 



One can also easily compute the thermal correction to the term TrM + , just by evaluating 
a tadpole of (j) fields, taking into account the SU{N) relation {A a , A 6 } = ^5 ab + 2d abc \ c . 
Finally, one can write the complete one-loop thermal effective action for soft modes as 

Z 2 = S 2 + 6S 2 , T = J d A x (~\f, u F^ - Y a ^ A ^) ( 2 - 25 ) 

- — I d A xd A yTr { r„ x (x) < x\ Q " Q \A y > r u \y)) 

12 J 4vr J y y -(Q ■ d) 2iy KyJ J 

- J d A x (fl(T) Tr{A%x)) + 32 Tr (h%(x) - H 2 (x) 



2 / Ar2 _ i rr2\ ^ 
Ti(M + (x)) 



fl L N 2 -IT 
4 ^ UN f 2 



where is a light-like four vector (see Ref. ||13||)) and |T5|-|T7 

24/1 



/ NT 2 \ 

/ir (T) = /.!--— , (2.26) 



and 

= M r„ - + [f „, r„] , (2.27) 

/ jVT 2 \ T 2 f 2 
<7(T) = C (l - --^J + -A . (2.28) 

If we put the external sources w M = a M = 0, and = Ql = Q, the third term in Eq. 
( |2.25| ) is just the HTL effective action ||18|| , which includes the Debye mass for the electric 
field. 

With the above effective action one can read the thermal correction to the quark con- 
densate, just by performing a functional derivative of Z 2 with respect to the external source 
s(x). It is given by 

/ 7V 2 -1T 2 \ 
(m) T = (qq) T=0 [l - -^Tp) > ( 2 - 29 ) 



in agreement with Ref. []15 



The value of the masses of the charged Goldstone bosons are obtained by finding the 
poles of the two-point functions of the axial current (see Ref. |§ for details about the 
computation). So one gets 
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(2.30) 



In the low T limit, and using the values of C(T) and f n (T) above one then finds, at 
leading order, 



(T) 




(2.31) 



If one puts C — 0, then the thermal masses of the charged bosons agree with those 



computed in Ref. |19| for scalar QED. Temperature corrections to the electromagnetic mass 
difference of pions have also been considered in We agree with the results of [20|, except 
for the fact that thermal corrections of order T 2 / f 2 were not considered there. 

Finally, let us stress that in the massless quark limit Dashen's theorem remains valid at 
finite T, as expected. 



III. SUM RULE APPROACH TO THE ELECTROMAGNETIC MASS 

DIFFERENCE OF PIONS 

A. Computation at Zero Temperature 

In this subsection we review the sum rule computation of the electromagnetic mass 
difference of pions. We will follow closely the notation and conventions of Ref. ||, working 
in this subsection in Minkowski space-time. In the following subsection we will generalize 
the computation at low T. 

At T = one works with the assumption that SU(2)r x SU(2)l is a symmetry of the 
QCD Hamiltonian which is spontaneously broken to SU(2)l + r. The pions are then the 
associated Goldstone bosons. 

In the chiral limit the neutral pion is massless, M^o = 0, while the charged pions get a 
mass from one-photon exchange contributions. The mass of the charged pions is given by 

M ^ = Tl d ' xA " V ^ {<P)\T{J»{x)J»m \<P)) , (3-1) 

where A Miy is the photon propagator, and is the electromagnetic current. In momentum 
space and in Feynman gauge, Eq. (ft.l|) becomes 



, ^ ie 2 r d 4 q m , 9 , 



2 J (2vr) 4 q 2 
where 

T^{q\vQ) I dtxe** (n(p)\T(J^x)J u (0))\n(p)) . (3.3) 
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Eq. ( p.3|) can be separated into non-contact (NC) and contact (C) contributions ||. The 
last one refers to one which has both photons interacting at the same vertex. 

One can reduce the non-contact contribution to the Compton scattering amplitude to 
vacuum polarization functions This is done using the soft pion theorem 



lim (7r a (p)P 

P^O \ 



O 



a 



--{(3 [Ql,0] 



a 



(3.4) 



where a and (3 are arbitrary states, and Q\ is the axial charge, and also the current algebra 



(3.5) 



In Eq. (|3.5|) and A a „ are the vector current and axial vector current, respectively. Thus, 
in the soft pion limit, 



limT(^)(g 2 ,p.g)^|(n^(g)-n^ 



|| /d 4 xe^(0 

J IT 



(3.6) 



T (VJ(*R 3 (0) - Al(x)Al(0) 



The above correlators are written in terms of spectral functions as 



n^ 3 (g) = iq 2 Wv- 



q 2 
q 2 



°° ds PV{S 







q 2 — s 



o 



00 , Pa(s) . f2 q^qy 
ds— if*-^- 



q 2 — s 



q- 



The contact term contribution to the Compton scattering amplitude is 

T^\q\p-q) = 2g^ . 
Adding the non-contact and contact contributions one gets 

Urn T {q 2 ,p- q) = - (g -StJ. 

p^o ^ \ q 



q" 



q 



(3.7) 
(3.8) 

(3.9) 

(3.10) 



Therefore, in the mass formula one has 

d A q 1 3ie 2 f d A q 



3ie z 



(2tt)V f 2 J (2^ Jo 



q 2 — s 



(3.11) 



The above integrals are evaluated using dimensional regularization. Thus, the first term in 
Eq. ( |3.11| ) vanishes. The second one is logarithmically divergent, but the coefficient of the 
divergent piece vanishes due to the second Weinberg sum rule. If one chooses a different UV 
regulator, there are also quadratic divergences which cancel due to the first Weinberg sum 
rule [§. 

The second term in Eq. ( |3.11| ) can be easily evaluated if at this point one assumes that 
the spectral functions py and pa are dominated by the vector and axial vector mesons 
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p v (s) = F 2 p 5 (s - M; 
p A {s)=F 2 A 5(s-M] 



Using the relations 



F 2 Ml 



F 2 A M 2 A , 



F 2 = F 2 p - f 2 , 



which are derived from Weinberg's sum rules [|| 

ds ( Pv (s) - pa(s)) = f 2 , 



dss(p v (s) - Pa(s)) = 



one finds 



3<= 2 F 2 M 2 
J^M 2 -^ In— p - 
lQn 2 ^ f 2 ln M 2 A ' 



(3.12) 
(3.13) 

(3.14) 

(3.15) 
(3.16) 

(3.17) 



which is the classical result of Das et al. p|. 

Better estimates of the correlator of vector and axial vector currents, not necessarily 
assuming a narrow- width spectral function as in Eqs. ( |3.12| )-( |3.13| ), have been considered in 
the literature M, ETJ. 



B. Computation at Low Temperature 



In this subsection we generalize the previous sum rule computation at low T, working in 
Euclidean space-time and in the ITF. At low T one can assume that the main hypothesis that 
hold at T = are still valid, and that one only needs to consider small thermal corrections 
to the previous computation. 

At low T chiral symmetry is still spontaneously broken, although at high enough T it 
is supposed to be restored. We make the assumption that at low T one can use the same 
kind of soft pion theorems than at T — 0. In the order of the computation we are working 
knowledge of the T = Weinberg's sum rules will be sufficient for us. A generalization of 



Weinberg's sum rules at finite T can be found in Ref. |22 
We start from Eq. 



Aty, but now evaluated in a thermal bath at equilibrium. We will 
follow the same steps as in Sect. [Ill A| to reduce the thermal amplitude 



7r(p)r(J M (x)J„(0)) 7r{p) 



T 



(3.18) 



to a thermal polarization function. Thus, we use the low T generalization of the soft pion 
theorem Eq. (|3.4|) at leading order in the thermal correction, i.e., 



lim 



TT a {p){3 



o 



a 



r 



A CO 



P [Qt, o] 



a 



(3.19) 
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and the current algebra. It is important to realize that Eq. ( |3.19|) only holds at low T and at 
leading order in T 2 / f 2 . At order T 4 //. 4 there are two distinct pion decay constants |17[], due 



to the loss of Lorentz invariance, and therefore at that order Eq. ( |3.19| ) would be modified. 

The thermal amplitude Eq. (|3.18| ) is then written in terms of the thermal correlators of 
vector and axial vector currents. In Euclidean space-time, these correlators are given by 

Tl v ^(K,T) = [ d i xe lK - x Y.{n\TV;(x)V a Me^' T \n) , (3.20) 

U*> a (K,T) = f d 4 xe^Y,(^^(x)A a u (0)e^/ T \n) , (3.21) 

J n 

where the sum is over the full set of eigenstates of the Hamiltonian H, and e~ n = 
J2 n (n\e~ H / T \Ti), and K is the Euclidean momentum. These thermal correlators have been 
computed in Ref. [T2[] assuming that the lightest particles of the spectrum of H, that is, 
the pions, give the main contribution to the sum. Then, just by using the same soft pion 
theorems that hold at T = and the current algebra, and after integrating over the pionic 
thermal space, one gets for massless pions 



U V (K,T) = [1-—, }U v (K,T = 0) + — rri(#,T = 0) , (3.22) 



T 2 \ y T 2 

-/"- ' \ ' gj2 J - T — 0) + jrj 2 "-\>"-- 

T 2 \ T 2 



U^(K,T) ={l-—j U^(K,T = 0) + — U^(K,T = 0) . (3.23) 

The two above low T correlators can thus be expressed in terms of their T = values, 
although there is a mixing between the vector and axial vector ones. Let us stress that at 
finite T, and due to the loss of Lorentz invariance in the thermal bath, the general form of the 
thermal correlator of two currents should depend on more unknown functions than the ones 
at T = 0. However, at low T one can assume that these correlators retain, approximately, 



their covariant form 12 



Once we know the value of the above thermal correlators, we are ready to compute the 
thermal corrections to the mass of the charged pions. We work in the ITF, where the thermal 
momentum measure is the one written in Eq. (|2.19| ). From Eqs. (|3.22|) - (|3.23|) we know that 




Ul u (K,T)-U^(K,T)=[l- — )(nl(K,T = 0)-UUK,T = 0)) . (3.24) 



In Feynman gauge, one then has 

M 2 ± (T)-3e 2 r**L±-J*-( 1 -H) [**L rdc Pv{s) - pA{s) (3 25) 

The thermal contribution of the first term in Eq. ( |3.25| ) does not vanish, as it happens for 
the T = contribution in dimensional regularization. Notice that we have neglected the 
thermal factors in front of the first term of Eq. ( |3.25| ) . Those terms would yield thermal 



corrections of order e 2 T 4 // 4 , which are subleading at the order of the computation that we 
are working. If one uses dimensional regularization to cure the UV divergencies of the above 
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integral, then the quadratic divergencies of Eq. ( |3.25|) vanish. If one uses a different UV 



regulator then to check that the quadratic divergencies which multiply the T 2 corrections 
cancel one would need the thermal factors that we have neglected in Eq. ( |3.25| ) . To simplify 
our computation, we will work with dimensional regularization, where Eq. ( |3.25| ) holds 
exactly. 

At this stage one can apply the hypothesis of vector meson dominance, which is valid at 
T = 0. We should then compute the integrals 



d 4 K ( F 2 F\ 



{2ti) a \K 2 + M 2 K 2 + M 2 A 



(3.26) 



After performing the sum over Matsubara frequencies one then gets the T = contribution 
(already evaluated in the previous subsection), plus the pure thermal part. However, since 
the masses of the vector resonances are much bigger than the temperatures we are consid- 
ering, M p , Ma 3> T, the thermal corrections associated to those particles are suppressed as 

~exp(-M PiA /^)- 

Thus, collecting all the thermal corrections, we get 

Ml ± (T) = (l - Ml ± + ^ , (3.27) 

plus subleading corrections. Therefore, we have checked our formula ( |2.30| ) for the N = 2 
case. 

It is worthwhile emphasizing that since the thermal corrections to the second term of 



Eq. (|3.25|) factorize, our results are also valid if we use a better evaluation of the spectral 



functions than the ones given in Eqs. (p3.12)-(|3.13). 



C. Results 



At zero temperature, the difference of the squared pion masses in the chiral limit is given 
by the result of Das et al., Eq. ( |3.17| ). Using the relations ( |3.14| ) to eliminate the parameters 



of the axial vector meson, A, it can be written as 



3a F 2 F 2 
M 2 ± -M 2 = -M 2 In , i 3.28 ) 



where a is the electromagnetic coupling constant. Although this result was obtained in the 
chiral limit, it can be used to give an estimate of the electromagnetic mass difference of 
pions. The difference of the squared pion masses is 



M; ± - M; ~ 2M 7r AAf 7r . (3.29) 

Taking the physical value of the pion mass M n = 135 MeV, and the following values for the 
other parameters involved: F w = 92.4 MeV, M p = 770 MeV, F p = 153 MeV and a = 1/137, 
one gets AM n = 4.8 MeV, which is in very good agreement with the experimental value 
{AM^^p = 4.6 MeV. 



12 



We thus expect that M 2 ± (T) — M 2 (T) at low temperature is also well approximated by 
the chiral limit calculation Eq. ( ft.27|) , i.e., 



AM^(T) = Ml ± {T) - Mlo{T) = 1 - — AM^ + vraT 2 , (3.30) 




where AM 2 = M 2 ± - M 2 o, given in Eq. ( |3.28|) . The result is plotted in Fig. 1, as a function 



of the temperature. Although our computation is only strictly valid for low T, we have 
extrapolated it up to T ~ V6f n ~ 220 MeV. The dashed-dotted line in Fig. 1 represents 
the part proportional to the electromagnetic mass difference at T = 0, AM 2 , in Eq. ( |3.30| ). 
and we see that it decreases with T. The dashed line corresponds to the typical thermal 
mass of charged particles, which grows with T, while the solid line is the full result. 



IV. CONCLUSIONS 

The difference of correlators of vector and axial vector currents can be taken as an order 
parameter of chiral symmetry breaking [23|,24]. The sum rule 



fl s ab = y d 4 x (o|r (v;(x)v» - ^(x)^(o) 



implies the asymmetry of the vacuum, since the operator V^(x)V MCt (0) — A^(x)A fia (0) trans- 
forms as the irreducible representation (3,3) of the symmetry group SU(2)l x SU(2)r, 
and therefore if the vacuum were invariant, the above correlator should be zero. If chiral 
symmetry is restored at high T, this order parameter should then vanish. 

Low T computations of the vector and axial vector correlators show a tendency towards 



chiral symmetry restoration |[L2]j . Although those computations cannot be trusted at higher 
T, they give a good estimate of the critical temperature of the phase transition, T c ~ 160 
MeV, which agrees with the values found in lattice computations . 

The electromagnetic mass difference of pions depends on the correlators of vector and 
axial vector currents 0. In concordance with the signs of chiral symmetry restoration at 
finite T, one would naturally expect that this mass difference should decrease with T. This 
is what we have obtained for low T, although it is not a priori obvious from our Eq. ( |3.30| ). 



We have found a slow decrease of this mass difference with T but only after a numerical 
analysis involving the T = values of the different parameters. It should be stressed that 
even if chiral symmetry is restored at a certain T c , typical thermal effects occurring in the 
plasma of pions and photons will always generate a thermal mass for the charged particles 
going as ~ eT. Therefore, one can never expect that the electromagnetic mass difference of 
pions should vanish at any high T. 

Our computation was done in the limit of massless quarks. In order to get a more 
trustable value of the mass difference (Mj^ — M|- )(T), one would need to take into account 
also QCD corrections proportional to the quark masses, apart from the pure electromagnetic 
correction computed here. This was, however, beyond the scope of this article. 
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APPENDIX A: HARD THERMAL LOOPS 

In a gauge field theory at finite T hard thermal loops are one-loop diagrams which are 
as important as the tree amplitudes for soft momenta. Soft denotes a energy scale <C T. If 
g is the gauge coupling constant, and it is assumed that g <C 1, then a soft scale is typically 
of the order ~ gT. For soft fields, then HTL's have to be resummed in order to take into 
account those one-loop effects consistently [[H]]. 

Braaten and Pisarski realized that HTL's only arise when all the external momenta of a 
diagram are soft and the internal one is hard ~ T. They were then able to establish some 
power counting rules to extract from each Feynman diagram the corresponding HTL. 



In the non-linear sigma model HTL's also appear for soft external momenta |14] , p!3| In 
%PT in the chiral limit without gauge interactions the expansion parameter for thermal 
corrections is T 2 /f 2 and it is this scale the one which allows to define the scales soft ~ 
T 2 / f 2 T or hard T. Resummation is not required in this case. 
In the chiral Lagrangian with electromagnetic interactions included there are then essen- 



tially two expansion parameters, e and yT 2 / f 2 . We have assumed in this article that the 
two of them are of the same order. 



In our case, one can then apply the same power counting rules derived in Ref. [[10 
Following Braaten and Pisarski's analysis for the Lagrangian ( |2.17|) , it is then easy to realize 
that vertices which do not carry a momentum dependence do not produce a HTL, except for 
tadpole-like diagrams. So apart from the tadpoles, the only vertices which produce HTL's 



arc 



■\F'A]) , (Ai) 

and 

- f -^Tt{d^H L )a^ . (A2) 
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FIG. 1. Electromagnetic mass difference of pions AM^(T) as a function of the temperature in 
the chiral limit. The dashed-dotted line corresponds to the piece proportional to the electromag- 
netic mass difference at T = 0, AM^, the dashed line to the term proportional to a which grows 
with T, and the solid line is the full result 
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